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Abstract 

A model of sound generated in a high subsonic 
(Mach 0.9) circular jet is solved numerically in 
cylindrical coordinates for nonaxisymmetric distur- 
bances. The jet is excited by transient mass injec- 
tion by a finite duration pulse via a rotating ring 
source. The flow field, near field and far field pres- 
sure disturbances corresponding to these sources are 
described. In particular, the resulting pressure field, 
which would serve to excite nearby panels, is illus- 
trated together with preliminary results on the ex- 
citation of thin slices of nearby panels. 

We consider both the short time behavior of the 
jet and the long time behavior, after the initial ex- 
citation pulse has exited the computational domain. 
The long time behavior of the jet is dominated by 
vorticity and pressure disturbances generated at the 
nozzle lip and growing as they convect downstream 
in the jet. These disturbances generate sound as 
they propagate. We find that rotating nonaxisym- 
metric disturbances persist for long times. Further- 
more, depending on location, both in phase and 
out of phase behavior can be found upon reflection 
across the jet axis. 

1 Introduction 

In this paper, the generation and propagation 
of sound in a high subsonic jet excited by nonax- 
isymmetric disturbances is simulated. We consider 
a circular jet exiting from a straight nozzle extend- 
ing to infinity in the upstream direction. The jet is 
excited by a nonaxisymmetric transient source and 
both the short time and long time behavior of the 
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disturbance pressure is considered. We solve the full 
nonlinear Euler equations in conservation form mod- 
ified so that a spreading jet profile, obtained from 
experimental measurements, is a solution in the ab- 
sence of any disturbances. 

The primary objectives are (i) to characterize non- 
axisymmetric aspects of the flow field and near field 
and (ii) to characterize nonaxisymmetric effects in 
the far field. The jet is initially excited by a spa- 
tially and temporally localized source of transient 
mass injection. We consider a rotating modulated 
ring source for the excitation, thus generating a ro- 
tating nonaxisymmetric disturbance which propa- 
gates along the jet in a helical fashion. This leads 
to an initial acoustic disturbance which propagates 
through the jet. As a result of the excitation, insta- 
bility waves are generated in the jet. These waves 
grow and then decay as they convect downstream, 
generating sound in the process. This phenomenon 
occurs over time scales much longer than that of 
the excitation pulse. The emphasis in this paper is 
on both the short time and the long time response 
of the jet, i.e., after the wave due to the excitation 
pulse has exited the domain of interest. Thus, insta- 
bility wave generated sound in the jet is simulated. 
The simulation does not directly account for sound 
generated by small scale turbulent sources in the jet. 

In previous work we have considered jet acous- 
tics for circular jets subject to both axisymmet- 
ric and nonaxisymmetric disturbances. In partic- 
ular, both linear and nonlinear responses have been 
considered. 1, 2> 3> 4 In these references the jet is ex- 
cited by a source of transient mass injection and 
the development of instabilities together with asso- 
ciated sound generation is computed. The results 
in these references demonstrate the properties of in- 
stability wave generated sound. A succession of dis- 
turbances is generated at the nozzle lip. These dis- 
turbances initially grow as they convect in the jet. 
As they grow they generate pressure disturbances 


1 





which propagate into the far field as sound. The 
far field radiation peaks at about 30° from the jet 
axis in agreement with measurements and observa- 
tions. Flapping modes, obtained from static modu- 
lated ring sources have also been computed. In these 
modes regions of large pressure alternate on different 
sides of the jet. 4 

We have also considered two dimensional Carte- 
sian jets. 5, 6 > 7 » 8 In two dimensions we have com- 
puted both the sound generated from instability 
waves and the interaction of this sound with a nearby 
flexible structure. These computations spanned a 
range of jet velocities up to and including super- 
sonic jets. The two dimensional results for high sub- 
sonic jets (the parameter range considered here) in- 
dicated that the long time behavior of the excited jet 
is dominated by a nearly periodic acoustic response 
at a particular fundamental frequency together with 
harmonics. 7 Peak radiation is beamed at approxi- 
mately 30° from the jet axis. The unsteady flow 
field is dominated by disturbances generated at the 
nozzle lip and propagating along the jet. In con- 
trast to the nearly periodic far field, the flow field 
exhibited an aperiodic, apparently chaotic behavior. 

The exact sources of jet noise have been 
identified from the basic equations of fluid 
dynamics. 9, 10, n * 12, 13 Computation of the exact 
sources would require solution of the Navier-Stokes 
equations in the jet accompanied by an appropriate 
turbulence model. Such computations are difficult 
to conduct, particularly in three dimensions. Both 
small scale and large scale disturbances in the jet 
can act as sources of sound. In this paper we con- 
sider only sound generated from large scale insta- 
bility waves. These instability waves are generated 
within the transiently excited jet by the mechanism 
described above. 

Instability waves or large scale structures which 
can act as sources of sound 
in a jet have been observed in experiments 1, 14, 15 
and studied by analytical 16, 17, 18, 19, 20, 21, 22 and 
numerical 1, ly 3 » 23 methods. While there have been 
a limited number of three dimensional computa- 
tions 24 ’ 25 , the cost of the 3D computations has to 
date precluded a detailed study of the role of non- 
axisymmetric large scale structures as sources of jet 
noise and has generally forced restrictive assump- 
tions such as linearity. On the other hand, many 
of the large scale structures observed in experiments 
are intrinsically nonaxisymmetric, e.g., helical and 
flapping modes, and thus a full understanding of jet 
acoustics requires an understanding of nonaxisym- 
metric modes in the jet and the resulting acoustic ra- 
diation. This is particularly true when the unsteady 


pressure in the jet serves to load nearby panels. In 
this case, the panel response can crucially depend on 
the nature of the nonaxisymmetric disturbances. 

In our model we solve a modified version of the full 
nonlinear Euler equations in cylindrical geometry. 
The equations are modified to account for a spread- 
ing jet as described below. As a result, the invis- 
cid, instability wave . sources of sound are computed 
directly together with the resulting sound genera- 
tion. We specifically include a straight nozzle from 
which the jet exits. The flow field and near field dis- 
turbances are computed together with the far field 
sound radiation. 

2 Model and Numerical Method 

The computational domain is shown in Fig- 
ures la. We consider a jet exiting from a straight, 
circular pipe of diameter D, where D — 4 inches. 
All boundaries are artificial boundaries at which ra- 
diation boundary conditions are employed. 1, 28 In 
these computations we employ an adaptive radia- 
tion boundary condition in that the assumed origin 
of coordinates is adaptively chosen so as to minimize 
a measure of the error of the boundary condition. 
The jet is excited by a transient source of mass in- 
jection located approximately ID downstream from 
the nozzle exit. This source is both spatially and 
temporally localized, leading to a broad band ex- 
citation. Unsteady pressure, density and velocity 
are computed both interior and exterior to the pipe. 
This leads to the generation of a train of instabil- 
ity waves which propagate along the jet, decaying 
beyond the potential core of the jet and generating 
sound. 

The primary objective is to consider the loading 
on nearby panels due to the jet pressure. Some lim- 
ited results along this line are presented below. In 
Figure lb we sketch the desired configuration. 

The Euler computations employ a (2-4) version of 
the MacCormack scheme. 26 Further details on the 
numerical scheme are given in the references. 5, 6j 7 
The spatial fourth order accuracy is necessary in or- 
der to follow the acoustic waves over large distances 
without undue numerical dispersion. This is partic- 
ularly important given the coarse grids necessitated 
by nonaxisymmetric computations. 

The Euler equations are solved in conservation 
form in cylindrical coordinates z,r,^>, where z de- 
notes the axial distance, r the radial distance and 
ip is the azimuthal angle. The dependent variable is 
the vector 

q= ( P,pu,pv,pw,E ) T , 
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where p is the density, u , v and w are the axial, radial 
and azimuthal velocity components respectively, and 
E is the total energy per unit volume. The pressure, 
p, is obtained from the equation of state. 

The Euler equations are modified in the jet do- 
main to account for the jet flow. The jet exits from 
a nozzle of diameter D located at z = 0 and the so- 
lution is computed both within and exterior to the 
nozzle. The Euler equations are modified to account 
for two different non-homogeneous forcing terms. 7 
One term serves as an excitation pulse to excite the 
jet. It corresponds to a localized source of mass in- 
jection. In axisymmetric computations we take a 
point source centered on the axis at z s . In the non- 
axisymmetric computations we consider a rotating 
ring source of radius r*, where r 8 ~ .3D. For each 
0 the source is a peaked Gaussian centered at r a , z 3 . 
The amplitude of the source is modulated The tem- 
poral dependence of the source is described in the 
references, e.g., 5 . The amplitude of the source is 
modulated by the function 1 + ei cos(0 + u>t) where 
ei — 0.125. Thus, the modulation represents a spin- 
ning excitation and serves to excite a spinning mode 
in the jet. The source is localized in time, thus lead- 
ing to broad band excitation. 

The second forcing term is designed so that in the 
absence of the excitation pulse the solution to the 
Euler equations would be a stationary profile cor- 
responding to a spreading jet. In the present com- 
putation, mean density and pressure are assumed 
constant, as axe the mean radial and azimuthal ve- 
locities. The mean axial velocity U is a taken from 
experimental measurements and models an isother- 
mal jet with exit Mach number 0.9. 27 

The inclusion of this term separates the compu- 
tation of the disturbance, in particular the result- 
ing instability waves, from the computation of the 
mean flow (i.e. the spreading jet). Thus, the result- 
ing system of equations allows for the simulation of 
instability waves and the resulting sound generation, 
together with the propagation of acoustic waves in 
the jet flow field, without requiring the computation 
of the spreading jet itself. Although this is a simpli- 
fied model, it captures many of the observed features 
of instability wave generated jet sound and permits 
high resolution computation of the coupling of jet 
noise with the flexible panels and the resulting radi- 
ation from the panels. In particular, the model al- 
lows for computation of many of the natural sources 
of jet noise (the instability waves) together with the 
sound radiated by these sources. 


3 Numerical Results 

We first consider the nature of the unsteady 
pressure in the excited jet. We describe the far field, 
near field and flow field disturbance pressure p = 
p(t, r, z, 0) - poo where poo is the ambient pressure. 
We do not consider a fully coupled response of the 
structure, so the unsteady pressure field, computed 
without the structure, will serve as the excitation 
pressure for the nearby panel. 

In response to the excitation source an initial 
acoustic wave emanates from the source location 
through the jet. Over short times this wave is 
the predominant feature in the unsteady flow field. 
However, this disturbance propagates essentially at 
the speed of sound and rapidly propagates out of the 
computational domain. A large slowly moving struc- 
ture, an instability wave, is generated by the excita- 
tion pulse and propagates downstream along the jet 
axis. The long time behavior of the jet is dominated 
by a shedding of disturbances from the nozzle lip, 
the propagation of these disturbances downstream 
and the generation of sound from these disturbances 
as they grow and decay along the jet. We note that 
there are also disturbances that propagate upstream 
inside the pipe. 

We first illustrate the structure of the fluctuating 
pressure field. Figures 2 and 3 are color contour plots 
for p in the r — z plane sin(0) = 0, i.e., the pressure 
field for the angles 0 = 0 and 0 = n. The coloring is 
according to the spectrum with red indicating large 
positive values of p and violet indicating large nega- 
tive values of p. Note that colors are assigned based 
on a logarithmic scale. These contours are shown in 
Figures 2 and 3 at times t\ = 14.8 and £2 = 10*6 
respectively, where the nondimensional time t is de- 
fined by t = tCoo/L , where t is dimensional time, Cqq 
is the ambient sound speed and the reference length 
L is 1 foot (3D). 

In interpreting these figures, we recall that the 
nozzle exit is at z — 0 and the initial source was a 
ring source centered approximately ID downstream 
of the nozzle exit and of approximately 0.5D in ra- 
dius. 

The initial acoustic disturbance due to the exci- 
tation pulse is still present in Figure 2. In addition, 
there is an array of acoustic disturbances propagat- 
ing upstream both above and below the jet. Acous- 
tic disturbances propagating up the nozzle can also 
be seen. These disturbances travel slowly, due to 
the nearly sonic velocity (Mach number 0.9) within 
the pipe. There is an array of large scale structures 
propagating along the jet axis with acoustic distur- 
bances emanating from them. These disturbances 
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grow and then decay as they propagate along the 
jet. 

At time £ 2 , corresponding to Figure 3, the initial 
acoustic wave has exited the computational domain. 
Additional large scale structures propagating along 
the jet axis have been generated, while those struc- 
tures in Figure 2 have propagated further down- 
stream, exhibiting a pronounced decay. In addition, 
there is a pronounced beaming of acoustic radiation 
at approximately 30° from the jet axis below the jet 
(i.e., in the half plane ip = n). 

We next illustrate different aspects of the un- 
steady pressure field. We first consider the far field 
pressure. In Figure 4 we plot p as a function of t for 
two different far field points. Both points are taken 
at 30° from the jet axis, however they are at opposite 
sides of the jet. Figure 5 is similar except that we 
consider points at 90° from the jet axis. Both figures 
are of the full time trace, including both the early 
and late time pressure. However, the predominant 
feature in the figures is the early time pressure cor- 
responding to the excitation source. The pressure is 
greater below the jet, i.e., at ip = n- This indicates 
that the predominant radiation from the excitation 
pulse occurs when the ring source has rotated so that 
maximum radiation occurs in the ip — 7r half plane. 

We next focus on the late time behavior of the jet 
by considering the time interval 12 < t < 20 after 
the initial pulse has exited the computational do- 
main. We plot in Figure 6 the time history over this 
interval. There is an irregular oscillation in the pres- 
sure. We note that this oscillation is out of phase 
between ip = 0 and ip = 7r. Furthermore, the differ- 
ence in level between ip — 0 and ip = 7T is no longer 
the same as for the initial pulse. Indeed, the overall 
sound pressure level, obtained from integrating p 2 
in time is higher at ip = 0 than for ip — n. This is 
shown in Figure 7 where the overall sound pressure 
level over this time interval is plotted against the far 
field angle for the half planes ip = 0 and ip = it. The 
data in the figure is normalized to 0 db for 90° for 
ip — 0. 

Thus the effect of the rotating source is to signif- 
icantly enhance the differences in the initial pulse, 
depending on what plane has a maximal source at 
the time of maximal radiation. However, nonax- 
isymmetric behavior in both amplitude and phase 
persists to late times after the excitation pulse has 
exited the computational domain. 

We next consider the near field pressure. In Figure 
8 we plot the late time pressure at approximately 3 D 
downstream and approximately 5 D from the jet axis. 
In contrast to the far field pressure, the pressure os- 
cillations are now in phase for the half planes ip — 0 


and ip = 7r. There is also a long wavelength modu- 
lation in the relationship between the amplitude of 
p between the two half planes so that initially p in 
the half plane ip = 7r is greater than at the corre- 
sponding point in the half plane ip — 0. The relative 
amplitudes then switch so that the pressure in the 
half plane ip = 0 is greater and then they switch 
back again. We believe that this is a residual effect 
of the rotation of the excitation source. 

We next consider the flow field pressure. In Figure 
9 we plot p at approximately six diameters down- 
stream and one diameter from the jet axis. There 
is a nearly sinusoidal oscillation which is again in 
phase for the two half planes. There is again a long 
time scale modulation for the relative amplitudes of 
the pressure between the 2 half planes. 

We next illustrate the persistence of the rotating 
nature of the disturbance in Figure 10 where we plot 
the axial velocity w as a function of the cylindrical 
angle ip for a range of values of t (12 < t < 14.5). 
The data is taken at approximately 1.5D down- 
stream at a distance of .5D from the jet axis. For 
these times the excitation pulse has long since de- 
cayed to zero. The figure shows an axial structure 
for w, which has completed approximately one full 
rotation within the time interval considered, thus 
illustrating the rotating nature of the long time so- 
lution. 

We next consider the effect of the unsteady pres- 
sure on nearby flexible panels. In previous two di- 
mensional computations we have considered a fully 
coupled model where both the panel response and 
radiation are computed simultaneously with the evo- 
lution of the unsteady pressure in the jet and are 
coupled to the computation of the pressure field in 
the jet. In view of the computational expense of the 
full three dimensional computations presented here, 
we consider a less complete model in which the pres- 
sure in the jet is first computed without the presence 
of the panel and is then used as a loading for a com- 
putation of panel response and radiation. Thus, the 
coupling between the panel radiation and the jet is 
neglected. 

The loading pressure is computed as part of the 
jet computation described above. At present, we 
have only considered slices of 8 inch panels located 
along the axis of the jet at a height of 3 D from the 
jet axis. We have computed the panel response to 
loading from the jet, employing parameters typical 
of aluminum. We illustrate this in Figure 11 where 
we plot the panel velocity up as a function of t for 
two panels centered at z = 0D and Z = 6D, re- 
spectively. There is a delayed response for the more 
downstream panel since it is at a great distance from 
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the jet source region. Furthermore, there is an in- 
creased emphasis of lower frequencies as the panel 
location shifts downstream. 

Finally, we illustrate in Figure 12 the rotation of 
pressure disturbances in the jet where we show cross 
sections of the unsteady pressure as a function of r 
and ^ at a fixed z location (approximately 2D down- 
stream from the nozzle). The diameter of the circu- 
lar cross sections is approximately ID. The cross 
sections are shown for four different times, t a = 7.5, 
U = 9.3, t c = 11.2, t d — 13.0, all of which are long 
after the excitation pulse has decayed to zero. The 
rotation of the pressure disturbance near the center 
of the jet can be clearly seen. We note that the rota- 
tion presented here is a different mode from the flap- 
ping modes described previously, where the nonax- 
isymmetry in the excitation source was stationary. 4 

4 Conclusion 

We have developed a code to solve the nonlinear 
Euler equations for nonaxisymmetric disturbances 
in an axisymmetric jet. The code is second order 
in time and fourth order in space and is written in 
cylindrical coordinates with the cylindrical axis the 
axis of the jet. The simulation includes both the 
near and far field of the jet as well as the upstream 
region including a straight pipe from which the jet 
exits. The jet is excited by a transient ring source 
located near the nozzle exit with a rotating, nonax- 
isymmetric modulation. The long time behavior of 
the jet is dominated by a train of disturbances gen- 
erated at the nozzle lip and propagating downstream 
along the jet. These disturbances in turn generate 
sound as they grow and decay in the jet shear layer. 

The rotating nature of the unsteady flow field per- 
sists to late times, after the initial excitation pulse 
has decayed to zero. We have compared the distur- 
bance pressure at different location in the computa- 
tional domain, focusing on differences between the 
half planes xp = 0 and xp = ir where xp is the cylin- 
drical angle. The differences between these two half 
planes for the long time pressure depends on loca- 
tion. In the far field there is almost a 180° phase 
difference for the long time pressure oscillations be- 
tween the two half planes. For near field locations, 
the pressure oscillations are in phase. For all points 
examined, there is a long time scale modulation in 
which the maximum pressure switches from one half 
plane to the other. 

Preliminary investigations of the use of the com- 
puted jet pressure as loading for an array of flexible 
panels are presented. These results demonstrate the 


increased low frequency response as the downstream 
location of the panels increases. 

Finally, examination of cross sections transverse 
to the jet axis indicates rotating pressure distur- 
bances. These rotating disturbances are different 
from the flapping disturbances we have previously 
computed using excitation sources where the nonax- 
isymmetry was stationary. 
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Figure 2. Color contour plot for p, t = 14.8 
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Figure 3. Color contour plot for p, t = 1 6.6 
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Figure 12. Axial cross sections of p at four successive times (times indicated in text.) 






